***************************************************************************
************* DESCRIPTIVE STATISTICS AND BALANCE TESTS ********************
***** Local Context and Mobilization in Poor Communities ******************
***************** Prisca Jöst and Ellen Lust ******************************
***************************************************************************

clear all

set maxvar 32000

use "REPLICATION_DATA_FINAL.dta", clear

*********************************************
**** PREPARATION *******************
*********************************************

*How many SQKM have n>20?
bysort sqkm: gen sqkm_freq= _N

*Which are those?
gen less20sqkm= 0 if sqkm_freq!=.
replace less20sqkm=1 if sqkm_freq<20

encode sqkm, gen(n_sqkm)

*Now exclude those from the data
drop if less20sqkm==1


**** Countries: Generate country variable ****

*** Zambia ****

gen Zambia_c = .
replace Zambia_c = 1 if LGPI_region=="Zambia" | LGPI_region=="Lusaka"
replace Zambia_c = 0 if LGPI_region=="Malawi" | LGPI_region=="Lilongwe" | LGPI_region=="Nairobi"

*** Malawi ****

gen Malawi_c = .
replace Malawi_c = 1 if LGPI_region=="Malawi" | LGPI_region=="Lilongwe"
replace Malawi_c = 0 if LGPI_region=="Zambia" | LGPI_region=="Lusaka" | LGPI_region=="Nairobi"

** Kenya ***
gen Kenya = .
replace Kenya = 1 if LGPI_region=="Nairobi"
replace Kenya = 0 if LGPI_region=="Zambia" | LGPI_region=="Lusaka" | LGPI_region=="Malawi" | LGPI_region=="Lilongwe"


*Social ties: Know some or most people in the neighborhood/village 

gen knowpeople= 0 if vnin_q3!=. & vnin_q3!=5
replace knowpeople= 1 if (vnin_q3==4|vnin_q3==3)

*Insufficient income 

gen insuffinc= 0 if wealth1!=. & wealth1!=5
replace insuffinc=1 if (wealth1==3 | wealth1==4)


******* SQKM MEASURES (CONTEXTUAL VARIABLES) ***********

*Here the measures on the sqkm level will be calculated: Share of Respondents who know some or most people 

bysort sqkm: egen knowpeopleshare= mean(knowpeople)
histogram knowpeopleshare, frequency

*Density: share by sqkm using 80, 60, 70 and 85 percent threshold

gen knowpeoplesqkm_binary= 0 if knowpeopleshare!=.
replace knowpeoplesqkm_binary= 1 if knowpeopleshare>0.79

gen knowpeoplesqkm_binary2= 0 if knowpeopleshare!=.
replace knowpeoplesqkm_binary2= 1 if knowpeopleshare>0.59

gen knowpeoplesqkm_binary3= 0 if knowpeopleshare!=.
replace knowpeoplesqkm_binary3= 1 if knowpeopleshare>0.69

gen knowpeoplesqkm_binary4= 0 if knowpeopleshare!=.
replace knowpeoplesqkm_binary4= 1 if knowpeopleshare>0.84


* Insufficient income by SQKM

bysort sqkm: egen sqkinsuffinc= mean(insuffinc)

histogram sqkinsuffinc, frequency

gen sqkminsuffinc_binary= 0 if sqkinsuffinc!=.
replace sqkminsuffinc_binary=1 if sqkinsuffinc>0.79

**************************************************
*************** INDIVIDUAL CONTROLS ***************
**************************************************

* CtrlAge

gen female = .
replace female = 1 if gender == 2
replace female = 0 if gender == 1
label define female 1 "female" 0 "male"
label values female female

** education **
gen educ2 = .
replace educ2 = 0 if education == 1
replace educ2 = 1 if education > 1
label define educ2 0 "no schooling" 1 "at least primary schooling"
label values educ2 educ2

** lived **
gen lived2 = .
replace lived2 = 1 if lived == 1
replace lived2 = 2 if lived == 2 | lived == 3 | lived == 4
replace lived2 = 3 if lived == 5
label define lived2 1 "less than 1 year" 2 "more than 1 year" 3 "all my life"

*** Asked by village head or neighbor *****

gen neighbor_VH = 0 if caul_q1<=5
replace neighbor_VH = 1 if caul_q1==1 | caul_q1==5

** Asked by village head ****

gen VH = 0 if caul_q1<=5
replace VH = 1 if caul_q1==1

*** Asked by neighbor ***
gen neighbor = 0 if caul_q1<=5
replace neighbor = 1 if caul_q1==5

*** Asked by MP ***
gen MP = 0 if caul_q1<=5
replace MP = 1 if caul_q1==4

*** Asked by TA ***
gen TA = 0 if caul_q1<=5
replace TA = 1 if caul_q1==2

*** Asked by councilor ***
gen councilor = 0 if caul_q1<=5
replace councilor = 1 if caul_q1==3

** encode regions variable ***

encode LGPI_region, gen(n_LGPI_region)



********************************************************************************
*Drop if insuficient income (poor) == 0 to only look at the poor in the sample*
********************************************************************************

drop if insuffinc == 0

*********************** DESCRIPTIVE STATISTICS: Appendix D ********************

***************************************************
** Motivations **
***************************************************

*** NOTE: RUN ANALYSIS BY COUTNRY AND FOR ALL POOR RESPONDENTS COMBINED
*drop if Kenya==0
*drop if Malawi==0
*drop if Zambia==0

*** Participation in Educational Programs and Funds (Figure D.1-3)

*Since {0} of last year have you or your household...Participated in construction/maintenance of school buildings or associated infrastructure (desks, chairs, etc.)
gen educ_contr = .
replace educ_contr = 1 if educ_q93==2
replace educ_contr = 0 if educ_q93==1

tab educ_contr

*Since {0} of last year have you or your household...Donated any school supplies
gen educ_contr2 = .
replace educ_contr2 = 1 if educ_q94==2
replace educ_contr2 = 0 if educ_q94==1

tab educ_contr2 


*Since {0} of last year have you or your household...Assisted with the school feeding program
gen educ_contr3 = .
replace educ_contr3 = 1 if educ_q95==2
replace educ_contr3 = 0 if educ_q95==1

tab educ_contr3 


*Since {0} of last year have you or your household...Helped provide security for the school
gen educ_contr4 = .
replace educ_contr4 = 1 if educ_q96==2
replace educ_contr4 = 0 if educ_q96==1

tab educ_contr4 

*Since {0} of last year have you or your household...Offered additional contributions
gen educ_contr5 = .
replace educ_contr5 = 1 if educ_q97==2
replace educ_contr5 = 0 if educ_q97==1

tab educ_contr5 

* Have you or your household made any of these contributions
gen any_contr = .
replace any_contr = 1 if educ_contr==1 | educ_contr2==1 |educ_contr3==1 |educ_contr4==1 |educ_contr5==1
replace any_contr = 0 if educ_contr==0 & educ_contr2==0 & educ_contr3==0 & educ_contr4==0 & educ_contr5==0

*** Reasons to Participate in Educational Fund ***

*Would you say that you participate,at least partly, because you can enjoy the event and company of others* educ_q106
gen educ_ben = .
replace educ_ben = 1 if educ_q106==2
replace educ_ben = 0 if educ_q106==1

tab educ_ben knowpeoplesqkm_binary, chi2
tab educ_ben knowpeoplesqkm_binary, col

*Would you say that you participate,at least partly, because you hear news and gain information* educ_q107
gen educ_ben2 = .
replace educ_ben2 = 1 if educ_q107==2
replace educ_ben2 = 0 if educ_q107==1

tab educ_ben2 knowpeoplesqkm_binary, chi2
tab educ_ben2 knowpeoplesqkm_binary, col

*Would you say that you participate,at least partly, because you or your household get goodies and gifts, or other material rewards* educ_q108
gen educ_ben3 = .
replace educ_ben3 = 1 if educ_q108==2
replace educ_ben3 = 0 if educ_q108==1

tab educ_ben3 knowpeoplesqkm_binary, chi2
tab educ_ben3 knowpeoplesqkm_binary, col


*Would you say that you participate,at least partly, because Others will think poorly of you or your household if you don't participate

gen educ_sanc = .
replace educ_sanc = 1 if educ_q109==2
replace educ_sanc = 0 if educ_q109==1

tab educ_sanc knowpeoplesqkm_binary, chi2
tab educ_sanc knowpeoplesqkm_binary, col


*Would you say that you participate,at least partly, because You or your household will have to pay fines, lose property or suffer other material loss if you don't
gen educ_sanc2 = .
replace educ_sanc2 = 1 if educ_q111==2
replace educ_sanc2 = 0 if educ_q111==1

tab educ_sanc2 knowpeoplesqkm_binary, chi2
tab educ_sanc2 knowpeoplesqkm_binary, col


*Would you say that you participate,at least partly, because You will be physically punished if you don't
gen educ_sanc3 = .
replace educ_sanc3 = 1 if educ_q113==2
replace educ_sanc3 = 0 if educ_q113==1

tab educ_sanc3 knowpeoplesqkm_binary, chi2
tab educ_sanc3 knowpeoplesqkm_binary, col


*Would you say that you participate,at least partly, because you personally think it is the right thing to do.
gen educ_sanc4 = .
replace educ_sanc4 = 1 if educ_q115==2
replace educ_sanc4 = 0 if educ_q115==1

tab educ_sanc4 knowpeoplesqkm_binary, chi2
tab educ_sanc4 knowpeoplesqkm_binary, col

*Thinking about the contribution of others, do you think that everyone contributes equally?
gen educ_equal = .
replace educ_equal = 1 if educ_q116==2
replace educ_equal = 0 if educ_q116==1

tab educ_equal knowpeoplesqkm_binary, chi2
tab educ_equal knowpeoplesqkm_binary, col


******** RECODE AS REWARDS, BENEFITS, RIGHT THING TO DO ********

** run analysis by country, drop other countries first **

*drop if Kenya==0
*drop if Malawi==0
*drop if Zambia==0

gen educ_rewards = .
replace educ_rewards = 1 if educ_ben==1 | educ_ben2==1| educ_ben3==1
replace educ_rewards = 0 if educ_ben==0 & educ_ben2==0 | educ_ben3==0

gen educ_sanctions = .
replace educ_sanctions = 1 if educ_sanc==1 | educ_sanc2==1| educ_sanc3==1 
replace educ_sanctions = 0 if educ_sanc==0 & educ_sanc2==0 & educ_sanc3==0 

gen educ_right = .
replace educ_right = 1 if educ_sanc4==1
replace educ_right = 0 if educ_sanc4==0


tab educ_rewards knowpeoplesqkm_binary, chi2
tab educ_sanctions knowpeoplesqkm_binary, chi2
tab educ_right knowpeoplesqkm_binary, chi2


*********************
**** Burial Fund *****
**********************

*Would you say that you participate, at least partly,because you can enjoy the event and company of others? 
gen educ_sanc21 = .
replace educ_sanc21 = 1 if extr_q25==2
replace educ_sanc21 = 0 if extr_q25==1

tab educ_sanc21 knowpeoplesqkm_binary if extr_q15_6==1, chi2
tab educ_sanc21 knowpeoplesqkm_binary if extr_q15_6==1, col


*Would you say that you participate, at least partly,because you hear news and gain information?
gen educ_sanc22 = .
replace educ_sanc22 = 1 if extr_q26==2
replace educ_sanc22 = 0 if extr_q26==1

tab educ_sanc22 knowpeoplesqkm_binary if extr_q15_6==1, chi2
tab educ_sanc22 knowpeoplesqkm_binary if extr_q15_6==1, col


*Would you say that you participate, at least partly,because you or your household get goodies and gifts, or other material rewards?
gen educ_sanc23 = .
replace educ_sanc23 = 1 if extr_q27==2
replace educ_sanc23 = 0 if extr_q27==1

tab educ_sanc23 knowpeoplesqkm_binary if extr_q15_6==1, chi2
tab educ_sanc23 knowpeoplesqkm_binary if extr_q15_6==1, col


*Would you say that you participate, at least partly,because others will think poorly of your or your household if you don't participate?
gen educ_sanc24 = .
replace educ_sanc24 = 1 if extr_q28==2
replace educ_sanc24 = 0 if extr_q28==1

tab educ_sanc24 knowpeoplesqkm_binary if extr_q15_6==1, chi2
tab educ_sanc24 knowpeoplesqkm_binary if extr_q15_6==1, col


*Would you say that you participate, at least partly,because you or your household will have to pay fines, lose property or suffer other material loss if you don't attend
gen educ_sanc25 = .
replace educ_sanc25 = 1 if extr_q30==2
replace educ_sanc25 = 0 if extr_q30==1

tab educ_sanc25 knowpeoplesqkm_binary if extr_q15_6==1, chi2
tab educ_sanc25 knowpeoplesqkm_binary if extr_q15_6==1, col


*Would you say that you participate, at least partly,because you will be physically punished if you don't
gen educ_sanc26 = .
replace educ_sanc26 = 1 if extr_q32==2
replace educ_sanc26 = 0 if extr_q32==1

tab educ_sanc26 knowpeoplesqkm_binary if extr_q15_6==1, chi2
tab educ_sanc26 knowpeoplesqkm_binary if extr_q15_6==1, col


*Would you say that you participate, at least partly,because you personally think it is the right thing to do.
gen educ_sanc27 = .
replace educ_sanc27 = 1 if extr_q34==2
replace educ_sanc27 = 0 if extr_q34==1

tab educ_sanc27 knowpeoplesqkm_binary if extr_q15_6==1, chi2
tab educ_sanc27 knowpeoplesqkm_binary if extr_q15_6==1, col


****** RECODE FOR REWARDS, SANCTIONS, RIGHT THING TO DO ****

** run analysis by country, drop other countries first **

*drop if Kenya==0
*drop if Malawi==0
*drop if Zambia==0

gen ex_rewards = .
replace ex_rewards = 1 if educ_sanc21==1 | educ_sanc22==1| educ_sanc23==1
replace ex_rewards = 0 if educ_sanc21==0 & educ_sanc22==0 & educ_sanc23==0

gen ex_sanctions = .
replace ex_sanctions = 1 if educ_sanc24==1 | educ_sanc25==1| educ_sanc26==1 
replace ex_sanctions = 0 if educ_sanc24==0 & educ_sanc25==0 | educ_sanc26==0 

gen ex_right = .
replace ex_right = 1 if educ_sanc27==1
replace ex_right = 0 if educ_sanc27==0


tab ex_rewards knowpeoplesqkm_binary, chi2
tab ex_sanctions knowpeoplesqkm_binary, chi2
tab ex_right knowpeoplesqkm_binary, chi2


***** Has your leder ever asked you to do this in the past? (Table D.1) *****

gen asked=.
replace asked=1 if caul_q18==2
replace asked=0 if caul_q18==1

*drop if Kenya==0
*** Voting ***
tab asked if caul_q2==1 & VH==1
tab asked if caul_q2==1 & neighbor==1
tab asked if caul_q2==1 & MP==1
tab asked if caul_q2==1 & TA==1
tab asked if caul_q2==1 & councilor==1

*** Burial Fund ***
tab asked if caul_q2==2 & VH==1
tab asked if caul_q2==2 & neighbor==1
tab asked if caul_q2==2 & MP==1
tab asked if caul_q2==2 & TA==1
tab asked if caul_q2==2 & councilor==1

*** Educational Fund ****
tab asked if caul_q2==3 & VH==1
tab asked if caul_q2==3 & neighbor==1
tab asked if caul_q2==3 & MP==1
tab asked if caul_q2==3 & TA==1
tab asked if caul_q2==3 & councilor==1

*drop if Malawi==0
*** Voting ***
tab asked if caul_q2==1 & VH==1
tab asked if caul_q2==1 & neighbor==1
tab asked if caul_q2==1 & MP==1
tab asked if caul_q2==1 & TA==1
tab asked if caul_q2==1 & councilor==1

*** Burial Fund ***
tab asked if caul_q2==2 & VH==1
tab asked if caul_q2==2 & neighbor==1
tab asked if caul_q2==2 & MP==1
tab asked if caul_q2==2 & TA==1
tab asked if caul_q2==2 & councilor==1

*** Educational Fund ***
tab asked if caul_q2==3 & VH==1
tab asked if caul_q2==3 & neighbor==1
tab asked if caul_q2==3 & MP==1
tab asked if caul_q2==3 & TA==1
tab asked if caul_q2==3 & councilor==1

*drop if Zambia==0
*** Voting ***
tab asked if caul_q2==1 & VH==1
tab asked if caul_q2==1 & neighbor==1
tab asked if caul_q2==1 & MP==1
tab asked if caul_q2==1 & TA==1
tab asked if caul_q2==1 & councilor==1

*** Burial Fund ***
tab asked if caul_q2==2 & VH==1
tab asked if caul_q2==2 & neighbor==1
tab asked if caul_q2==2 & MP==1
tab asked if caul_q2==2 & TA==1
tab asked if caul_q2==2 & councilor==1

*** Educational Fund ***
tab asked if caul_q2==3 & VH==1
tab asked if caul_q2==3 & neighbor==1
tab asked if caul_q2==3 & MP==1
tab asked if caul_q2==3 & TA==1
tab asked if caul_q2==3 & councilor==1

***** Communities by Density of Social Ties Table D.2) ****

tab educ2 knowpeoplesqkm_binary, col row

tab female knowpeoplesqkm_binary, col row

tab lived2 knowpeoplesqkm_binary, col row

tab highPopulation knowpeoplesqkm_binary, col row

tab insuffinc knowpeoplesqkm_binary, col row

tab CtrlAge knowpeoplesqkm_binary, col row

***** Number of communities and respondents (Table D.3) *****

tab knowpeoplesqkm_binary

tab knowpeoplesqkm_binary if Kenya==1
tab knowpeoplesqkm_binary if Zambia==1
tab knowpeoplesqkm_binary if Malawi==1

codebook n_sqkm if Kenya==1 & knowpeoplesqkm_binary==1
codebook n_sqkm if Zambia==1 & knowpeoplesqkm_binary==1
codebook n_sqkm if Malawi==1 & knowpeoplesqkm_binary==1

codebook n_sqkm if Kenya==1 & knowpeoplesqkm_binary==0
codebook n_sqkm if Zambia==1 & knowpeoplesqkm_binary==0
codebook n_sqkm if Malawi==1 & knowpeoplesqkm_binary==0

codebook n_sqkm if knowpeoplesqkm_binary==1
codebook n_sqkm if knowpeoplesqkm_binary==0


*********** Descriptives by Community ***********
* Subsamples by countries and rural/urban areas
*drop if Kenya==1 | Malawi_c==1
*drop if Kenya==1 | Zambia_c==1
*drop if Zambia_c==1 | Malawi_c==1
*drop if highPopulation==1
*drop if highPopulation==0


*keep sqkm knowpeople 

*collapse (mean) knowpeople, by(sqkm)

*sum knowpeople 

******************************************************************
*********** BALANCE TESTS APPENDIX D (Tables D.11-15)***********
******************************************************************

*********** Dropped from experment because no leader (Table D.4)***********
tab caul_q5
* drop no and do not know answers: 1,216 plus 685= 1901 (13.47%)

**** by leader type ****
gen dropped =.
replace dropped = 1 if caul_q5==1 | caul_q5==3
replace dropped = 0 if caul_q5==2

		  **  1 Village Head/Neighborhood Block Leader 
          **  2 Traditional Authority 
           ** 3 Local Councilor 
           ** 4 Member of Parliament 
          **  5 next door neighbor 
		  
** Treatment occurance by gender, age, urban/rural, social density


** Gender **
tab neighbor_VH female if dropped==0, row col
tab neighbor_VH female if dropped==0, chi2

tab PGtreat female if dropped==0, row col
tab PGtreat female if dropped==0, chi2

tab burialtreat female if dropped==0, row col
tab burialtreat female if dropped==0, chi2

tab votetreat female if dropped==0, row col
tab votetreat female if dropped==0, chi2

tab commonitortreat female if dropped==0, row col
tab commonitortreat female if dropped==0, chi2

tab leadermonitortreat female if dropped==0, row col
tab leadermonitortreat female if dropped==0, chi2

tab fewtreat female if dropped==0, row col
tab fewtreat female if dropped==0, chi2

tab alltreat female if dropped==0, row col
tab alltreat female if dropped==0, chi2

** Education **
tab neighbor_VH educ2 if dropped==0, row col
tab neighbor_VH educ2 if dropped==0, chi2

tab PGtreat educ2 if dropped==0, row col
tab PGtreat educ2 if dropped==0, chi2

tab burialtreat educ2 if dropped==0, row col
tab burialtreat educ2 if dropped==0, chi2

tab votetreat educ2 if dropped==0, row col
tab votetreat educ2 if dropped==0, chi2

tab commonitortreat educ2 if dropped==0, row col
tab commonitortreat educ2 if dropped==0, chi2

tab leadermonitortreat educ2 if dropped==0, row col
tab leadermonitortreat educ2 if dropped==0, chi2

tab fewtreat educ2 if dropped==0, row col
tab fewtreat educ2 if dropped==0, chi2

tab alltreat educ2 if dropped==0, row col
tab alltreat educ2 if dropped==0, chi2

** Age **
tab neighbor_VH CtrlAge if dropped==0, row col
tab neighbor_VH CtrlAge if dropped==0, chi2

tab PGtreat CtrlAge if dropped==0, row col
tab PGtreat CtrlAge if dropped==0, chi2

tab burialtreat CtrlAge if dropped==0, row col
tab burialtreat CtrlAge if dropped==0, chi2

tab votetreat CtrlAge if dropped==0, row col
tab votetreat CtrlAge if dropped==0, chi2

tab commonitortreat CtrlAge if dropped==0, row col
tab commonitortreat CtrlAge if dropped==0, chi2

tab leadermonitortreat CtrlAge if dropped==0, row col
tab leadermonitortreat CtrlAge if dropped==0, chi2

tab fewtreat CtrlAge if dropped==0, row col
tab fewtreat CtrlAge if dropped==0, chi2

tab alltreat CtrlAge if dropped==0, row col
tab alltreat CtrlAge if dropped==0, chi2

** Lived **
tab neighbor_VH lived2 if dropped==0, row col
tab neighbor_VH lived2 if dropped==0, chi2

tab PGtreat lived2 if dropped==0, row col
tab PGtreat lived2 if dropped==0, chi2

tab burialtreat lived2 if dropped==0, row col
tab burialtreat lived2 if dropped==0, chi2

tab votetreat lived2 if dropped==0, row col
tab votetreat lived2 if dropped==0, chi2

tab commonitortreat lived2 if dropped==0, row col
tab commonitortreat lived2 if dropped==0, chi2

tab leadermonitortreat lived2 if dropped==0, row col
tab leadermonitortreat lived2 if dropped==0, chi2

tab fewtreat lived2 if dropped==0, row col
tab fewtreat lived2 if dropped==0, chi2

tab alltreat lived2 if dropped==0, row col
tab alltreat lived2 if dropped==0, chi2

** Density (Tables D.5-10)**
drop if dropped == 1

ttest knowpeopleshare, by(neighbor_VH) 
ttest knowpeopleshare, by(PGtreat) 
ttest knowpeopleshare, by(burialtreat) 
ttest knowpeopleshare, by(votetreat) 
ttest knowpeopleshare, by(commonitortreat) 
ttest knowpeopleshare, by(leadermonitortreat) 
ttest knowpeopleshare, by(fewtreat) 
ttest knowpeopleshare, by(alltreat) 

******* Are those who were dropped significanlty different from those who were not ****  
tab dropped female, row col
tab dropped female, chi2

tab dropped educ2, row col
tab dropped educ2, chi2

tab dropped CtrlAge, row col
tab dropped CtrlAge, chi2

tab dropped lived2, row col
tab dropped lived2, chi2

******* By Leader Type *******
		   *1 Village Head/Neighborhood Block Leader 
           *2 Traditional Authority 
           *3 Local Councilor 
           *4 Member of Parliament 
           *5 next door neighbor 
		   
* Village Head
tab dropped female if caul_q1==1, row col
tab dropped female if caul_q1==1, chi2

tab dropped educ2 if caul_q1==1, row col
tab dropped educ2 if caul_q1==1, chi2

tab dropped CtrlAge if caul_q1==1, row col
tab dropped CtrlAge if caul_q1==1, chi2

tab dropped lived2 if caul_q1==1, row col
tab dropped lived2 if caul_q1==1, chi2


* TA
tab dropped female if caul_q1==2, row col
tab dropped female if caul_q1==2, chi2

tab dropped educ2 if caul_q1==2, row col
tab dropped educ2 if caul_q1==2, chi2

tab dropped CtrlAge if caul_q1==2, row col
tab dropped CtrlAge if caul_q1==2, chi2

tab dropped lived2 if caul_q1==2, row col
tab dropped lived2 if caul_q1==2, chi2

* LC
tab dropped female if caul_q1==3, row col
tab dropped female if caul_q1==3, chi2

tab dropped educ2 if caul_q1==3, row col
tab dropped educ2 if caul_q1==3, chi2

tab dropped CtrlAge if caul_q1==3, row col
tab dropped CtrlAge if caul_q1==3, chi2

tab dropped lived2 if caul_q1==3, row col
tab dropped lived2 if caul_q1==3, chi2

* MP
tab dropped female if caul_q1==4, row col
tab dropped female if caul_q1==4, chi2

tab dropped educ2 if caul_q1==4, row col
tab dropped educ2 if caul_q1==4, chi2

tab dropped CtrlAge if caul_q1==4, row col
tab dropped CtrlAge if caul_q1==4, chi2

tab dropped lived2 if caul_q1==4, row col
tab dropped lived2 if caul_q1==4, chi2

* Neighbor
tab dropped female if caul_q1==5, row col
tab dropped female if caul_q1==5, chi2

tab dropped educ2 if caul_q1==5, row col
tab dropped educ2 if caul_q1==5, chi2

tab dropped CtrlAge if caul_q1==5, row col
tab dropped CtrlAge if caul_q1==5, chi2

tab dropped lived2 if caul_q1==5, row col
tab dropped lived2 if caul_q1==5, chi2
